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ON \ Abstract. The possibility of detecting the gravitomagnetic clock effect using 

artificial Earth satellites provides the incentive to develop a more intuitive approach to 
^ . its derivation. We first consider two test electric charges moving on the same circular 

s ! | orbit but in opposite directions in orthogonal electric and magnetic fields and show that 

p~ • the particles take different times in describing a full orbit. The expression for the time 

difference is completely analogous to that of the general relativistic gravitomagnetic 

f^) ' clock effect in the weak-field and slow-motion approximation. The latter is obtained by 

considering the gravitomagnetic force as a small classical non-central perturbation of 
the main central Newtonian monopole force. A general expression for the clock effect 
is given for a spherical orbit with an arbitrary inclination angle. This formula differs 
from the result of the general relativistic calculations by terms of order c~ 4 . 
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1. Introduction 

The general relativistic clock effect, as worked out in |l], |2|, consists in the difference in 
the orbital periods of two clocks moving in opposite directions along a circular equatorial 
orbit around a central rotating mass. It is not an easy task to derive the general effect 
within the framework of Einstein's theory of gravitation; in fact, for orbits of arbitrary 
inclination to the equatorial plane, the clock effect has been derived only for the case 
that the orbiting bodies describe spherical orbits of constant "radius" 0. Moreover, 
the case of elliptical orbits has not yet been investigated. 

The purpose of this paper is to provide an alternative derivation of this effect 
for orbits with zero eccentricity, in order to make its comprehension more intuitive by 
stressing and elucidating the analogies and differences between the weak-field and slow- 
motion approximation of general relativity and electromagnetism. General expressions 
are given for arbitrary values of inclination as well. They are useful in view of the recent 
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efforts devoted to exploring the possibility of measuring the clock effect by means of 
artificial near-Earth satellites || f|, |], |6], |7j . 

The paper is organized as follows: In Section 2 we deal with an electromagnetic 
system consisting of two test charges orbiting in opposite directions, acted upon by a 
central electric field and by a weaker orthogonal magnetic field. The latter is treated 
perturbatively to first order. Radiative and 0(v 2 /c 2 ) effects in such a system are 
neglected in order to outline the main features common to the gravitational case that 
is treated in Section 3. In this work the gravitomagnetic force is viewed classically as 
a non-central small linear perturbation of the main central Newtonian gravitoelectric 
monopole force. Section 4 contains concluding remarks. 

2. The electromagnetic scenario 

Let us consider two identical point charges q of mass m orbiting a central spherically 
symmetric distribution of total charge Q of opposite sign, e.g. Q < and q > 0. We 
suppose that the two charges follow identical but opposite closed circular orbits and 
denote the speeds of the counterclockwise and clockwise moving charges by v + and t>_, 
respectively, i.e. we assume the counterclockwise direction to be positive. In cylindrical 
coordinates {p, cf>}, the equation of motion of the two charges reads 

v 2 
m^- = \q\ E, (1) 

P 

with E = \Q\/p 2 and both charges describe a complete orbit in the same time T^ -* 



rw = 2*£ = 27T /W 



mp3 (2) 



v± n ~ \]\q\E \\qQ\' 

where it should be noted that the mean motion (i.e. orbital frequency) n = J\q\ E/mp 
depends on the charge-to-mass ratio q/m. 

If we switch on a magnetic field B = Bi z orthogonal to the plane of motion, the 
two charges will experience an additional Lorentz force F^ = (q/c)v x B, which for 
q > will be antiparallel to the electric field for the counterclockwise moving charge 
and parallel to the electric field for the clockwise moving charge. Hence the equation of 
motion for q > becomes 

2 

m^± = qE =f -v±B, (3) 

p c ~ 

and therefore 



2mc J \ Am 2 c 2 / 

We assume that the magnetic field is weak; therefore, we can neglect the square of 
the Larmor frequency q 2 B 2 /4m 2 c 2 in comparison with the square of the orbital frequency 
n 2 in eq. (||) and obtain v± — (n =p qB/2mc)p, i.e. 

(f)-(-H)< < 5 > 
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or 

dt ± = ± - ** B ~±**(i±J?-). (6) 

By integrating eq. (§) from to 2ir for the counterclockwise orbit and 2ir to for the 
clockwise orbit, we find the orbital periods of the two charges 

and hence their difference after one revolution 

2ixB 
T+-T.= -- P . (8) 

By inspection of eq. (^|) we see that the magnetic correction to the orbital period 
is independent of the charge-to-mass ratio of the orbiting electric charges, in contrast 
to the unperturbed period T^°\ Indeed we could have started Section 2 with the less 
restrictive assumption that the two charges only have the same charge-to-mass ratio. 

It is interesting to observe that eq. (H) is an exact consequence of eq. @, i.e. 
the weak-field approximation is unnecessary for the validity of this relation. This 
circumstance has an analogue in the gravitational case discussed in Section 3; that 
is, eq. (|17D below turns out to be exact for circular orbits in the equatorial plane of the 
exterior Kerr spacetime. 

Let us now assume that the two charges are far away from the central charge and 
current distributions so that the magnetic field can be considered to be generated by a 
magnetic dipole m = — pi z of magnitude /i = IS/c, where S is the surface area of the 
loop enclosed by the current / and therefore 

p 6 p z 

which, upon inserting into eq. ([£) yield 

T + - T - = -isp (10) 

c \Q\ 

This time difference depends on both the sign of the charges and on the direction of 
the magnetic field. Upon exchanging the signs of Q and q, i.e. Q > and q < 0, the 
counterclockwise revolving charge will move faster while the clockwise moving charge will 
move slower. However, as expected by charge symmetry, eq. ( |T0D will be unaffected if 
the signs of the charges and of the magnetic field are reversed simultaneously. Further, 
we note that the radius of the orbit does not appear in eq. (0) and that this time 
difference can be interpreted as a consequence of the fact that the speed of light has a 
finite value; indeed, for c — > oo, T + — T_ — > 0. 

The main considerations of this section are related, via the Larmor theorem, to 
certain interesting phenomena in rotating frames of reference ||, ||. 

Orbits of charged particles off the equatorial plane are no longer spatially closed 
in general because their instantaneous planes undergo Larmor precession induced by 
the magnetic field and it will thus be necessary to define the relevant periods T± in 



An alternative derivation of the gravitomagnetic clock effect 4 

terms of azimuthal closure involving a complete loop in the coordinate. For the 
description of such a configuration it is useful to introduce the local frame attached 
to the moving particle, where ir, It, ijv denote its orthogonal unit vectors related to 
the radial, along-track and cross-track directions, respectively. In fact, ir denotes the 
orthogonal direction, in the instantaneous orbital plane, to the radial one; in general, 
it does not coincide with the direction along the track unless one considers circular 
orbits. We will return to the case of spherical orbits in Section 4, where the local frame 
described here is employed along the unperturbed orbit. 

3. The gravitational scenario 

The electromagnetic scenario previously described is analogous to the following 
gravitational one: Let us consider a central spherically symmetric mass M rotating 
slowly with its proper angular momentum directed along the z-axis, J = Ji z , of 
an asymptotically inertial frame K{x, y, z} whose (x, y)-plane coincides with the 
equatorial plane of the gravitating source, and a pair of test bodies orbiting along a 
circular equatorial path in opposite directions. Further, we assume that the radius p 
of the orbits is much larger than the Schwarzschild radius r g = 2GM/c 2 of the central 
body; this would apply to an experiment in the field of the Earth, for instance. It is 
well known that in the weak-field and slow-motion approximation of general relativity 
the stationary spacetime metric of a rotating "spherically symmetric" mass-energy 
distribution generates the so-called gravitoelectric and gravitomagnetic fields (JTTJ|, |TT | 



GM. , N 

E„=-— in, (11) 






1G 
B g = V x A g = — 

c 



J - 3(J ■ ifl)i fl 






(12) 



with the gravitomagnetic potential given by 

A. = -*?i4I. (13) 

c r 6 

At this point the calculations follow closely those of the electromagnetic case previously 
examined in Section 2 because the equation of motion of a test body in the "weak" 
gravitational field of a general stationary axisymmetric mass-energy distribution is 
analogous to that of a point charge q acted upon by electric and magnetic fields E 
and B, 

maj = m(E 9 + -xB 9 ). (14) 

Therefore, by reasoning as in the electromagnetic case for circular orbits with zero 
inclination (Section 2), we obtain 

T+ - r_ = *& (15) 

C Eg 
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For equatorial circular orbits, eq. (|T2| ) immediately yields 

B S = f £.. (16) 

By inserting i? s and B g from eqs. ( |TT]) and ( |16D into eq. ([H]) we obtain the well-known 
expression 

t -- t -=^4m- < i7 > 

It is an interesting feature of eq. ( |TTD that the mass moving in the same sense 
of rotation as the central mass moves slower than the mass moving in the opposite 
direction. If we reversed the sense of rotation of the central gravitating source, the 
clockwise moving test mass would be slower. In this way the sense of rotation of the 
central mass is no longer a matter of convention but could be related to a physical 
phenomenon, i.e. the mass loop moving slower. Also in this case, in the limit c — > oo, 
T+-T- -> 0. 

An interesting feature of gravitoelectromagnetism is the gravitational Larmor 
theorem [|l!| according to which gravitoelectric and gravitomagnetic fields are locally 
equivalent to translational and rotational accelerations of an observer in Minkowski 
spacetime. The gravitational Larmor theorem is therefore in effect Einstein's principle 
of equivalence in the gravitoelectromagnetic context. Note that in eq. (|17]) the 
gravitoelectric Keplerian periods T^ = 2n/n = JGM/p 3 cancel out, no matter what 
the masses of the orbiting particles are, in accordance with the equivalence principle 
and in contrast to the electromagnetic case, where the unperturbed periods depend on 
the charge-to-mass ratio of the particles. 

Finally, it is worth noting that the gravitomagnetic correction to the unperturbed 
Keplerian period is independent of the radius of the orbit, a feature also common to 
the electrodynamic case. Moreover, the Newtonian constant G does not appear in eq. 
(|17|): this fact may account for the unexpectedly large value of the time shift in the field 
of the Earth which amounts to about 10~ 7 s. In some sense, this classical effect is a 
gravitomagnetic analog of the topological Aharonov-Bohm effect. It turns out that if we 
deal with the proper periods of the test bodies, i.e. the periods according to comoving 
clocks, then the result is the same as in eq. (17) up to terms of order c -4 that depend 
on G and p. 

4. Arbitrary inclination: inertial azimuthal closure 

If the orbital plane has an arbitrary inclination i to the equatorial plane of K{x, y, z}, 
a similar reasoning as in Section 2 holds: the orbital plane undergoes Lense-Thirring 
precession, which is the gravitational analog of the magnetic Larmor precession. Let us 
imagine that in the absence of the gravitomagnetic field, the orbit is a circle of radius 
r in a fixed plane that is inclined with respect to the equatorial (x, ?/)-plane by the 
inclination angle i; that is, the normal to the orbital plane is tilted away from the z-axis 
by the angle i. Moreover, the longitude of the ascending node is given by the azimuthal 
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angle Q, as in Figure 1. Once the gravitomagnetic field is "turned on," the orbit will 
twist out of this fixed plane. 

In order to derive analytically the time T needed to pass from O at to = to 
0o + 27r at to + T in the equatorial plane of the inertial observer we will use the following 
reference frames: the asymptotically inertial frame K{x, y, z} previously defined and 
a frame K {X, Y, Z} with the Z-axis directed along the orbital angular momentum L 
of the unperturbed test body, where the (X, F)-plane coincides with the unperturbed 
orbital plane of the test particle and the X-axis is directed along the line of nodes. 
K{x, y, z} and K {X, Y, Z} have the same origin located at the center-of-mass of 
the central body. The transformation matrix H. x x for the change of coordinates from 
K {X, Y, Z} to K{x, y, z] is given by 

/ cos Q — sin Q cos i sin Q sin i 



R 



xX 



sinfi cosficosi — cosfisinz | . (1£ 

V sinz cosi 

Our calculation for the azimuthal period T will be valid for a counterclockwise orbit as 
in Figure 1; however, it is clear from the symmetry of the configuration that the result 
for the clockwise case (T_) can be simply obtained from our result by reversing the sign 
of the perturbation term in T+. In the asymptotically inertial frame K {X, Y, Z}, let 
us choose cylindrical coordinates {p, tp, Z} and write eq. ([14]) in these coordinates. 
The unperturbed orbit is given by p = r = constant, Z = and (p = (fo + nt (with 
n = JGM/p 3 ). Since we will consider the gravitomagnetic acceleration c _1 v x B 9 as 
a small perturbation of the main gravitoelectric monopole term E s , we evaluate the 
disturbing acceleration with respect to the unperturbed orbit. By using eqs. (|TT|) and 
(|12|), eq. (|H]) can be written as 



GM GJ 

7T + 2^-^ncosi, (19) 

pi c ipi 

a v = pp + 2pp = 0, (20) 

az = Z = —n 2 Z + 4-r— -nsinisin (<p + nt). (21) 

c 2 p 2 

Here the radial and along-track components of the acceleration are given by a p = 

Xcos (p + Ysin p and a^ = — Xsin p + Ycos <p, respectively. Note that in eq. ([2~l|) 

we retain Z, because due to the gravitomagnetic non-central acceleration, the motion is 

no longer confined to a plane and therefore Z will be proportional to the gravitomagnetic 

perturbation and treated to first order. 



For a spherical orbit p remains constant to first order; therefore, from eq. (19) we 

obtain 

GJ 

p> = cp + nt r-^-tcosi = p + Sip, (22) 

c z p- i 

where ip = ipo + nt. It follows from the linear dependence of ip on t in eq. (22) that eq. 

(20) is satisfied. The solution of eq. (|2~ID reads 

GJ 

Z = — 2— — t sin i cos (p>o + nt) + k\ sinnt + k 2 cos nt, (23) 

c 2 p 2 
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where k\ and k 2 are constants of integration. If we assume that at to — the perturbed 
orbit agrees with the unperturbed orbit, i.e. Z(0) = 0, then A; 2 = and we will treat 
k\ to first order of the perturbation in what follows. Let us write the general solution 
of eqs. (p^-pTI) as X = p cos <p(t) and Y = psmip(t) together with eq. (f23|). By using 
cos(p + 5p) ~ cos p — Sip sin p and sin(p + 5ip) ~ sin p + 5ip cos p, we find 

Ct J . . . . 

X = pcosp + -z-ztcosismp, (24) 

c l p l 

GJ 

Y = psinp — — tcosicosy?, (25) 

c 2 p 2 

GJ 

Z = — 2^r-^tsmicosp + k\ sin ni. (26) 

&pr 

From these equations we obtain by means of the transformation flIB|) the solution in 
K{x, y, z}, 



+ -^j£[cosicosfisin</> + (cos 2 i — 2sin 2 i) sin f2 cos <£] , (27) 



H — — — t[cosisinf2smy3 — (cos i — 2 sin i) cos fi cos p] , (28) 



z = p sin i sin p + k\ cos i sin nt — — - t sin 2i cosy. (29) 



x = p(cos Q cos ip — cos i sin O sin p) + fci sin % sin f2 sin nt 

GJ 

c 2 p 
y = p(sin O cos <£ + cos i cos fi sin </?) — k\ sin i cos Q sin nt 

GJ 
c 2 p 

3 GJ 
2^ 
The temporal behavior of the azimuthal angle can be obtained via 

tan0=^-|, (30) 

and the time T needed to pass from O at t = to 0o + 27r at to + T follows upon 
expanding the relation 

tan0 o = tan(0 o + 27r), (31) 

where 

, y(0) sin Q cos ipo + cos i cos f2 sin ip 

tan0 o = ^^ = p: • • n • ' ^ 

x(U) cos i 2 cos (^o — cos z sm i I sin ^o 

and 

tan(0 o + 27r) = 4^- (33) 

x(T) 

Since the deviation from the unperturbed Kepler period T^ = 2tt /n will be small, let 
us write 

T=-(l + e), (34) 

n 
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with e< 1 and further 

sin nT = sin(27r + 2ne) ~ 27re, (35) 

cos <p(T) = cos(<po + nT) = cos(<po + 2n + 2ne) ~ cos <po — 2ire sin <p ,(36) 
sin <p(T) = sin(<p + nT) = sin(<p + 2tt + 2ire) ~ sin <p + 2ne cos <p . (37) 

Therefore, in the calculations for x(T) and y(T), terms proportional to k\ sin nT, due 
to eq. fl35|), will be of second order and will be neglected. Hence we find 

x(T) = p(cos Q cos ifo — cos i sin Q sin ipo) — 2-ntp (cos Q sin <po + cos i sin fi cos <po) 

GJ 1 

+ 27T-5— - — [cos i cos f2 sin ip + (cos 2 z — 2 sin 2 z) sin fl cos ^ ] , (38) 

c z p 2 n 

y(T) = p(sin Q cos y?o + cos i cos f2 sin <^o) ~ 27cep (sin fi sin ipo — cos i cos Q cos y?o) 

GJ 1 
+ 2-k——- - — [cos i sin f2 sin o? — (cos 2 i — 2 sin 2 i) cos f2 cos o^ol • (39) 



Using eqs. fl3TD -(p3D and eqs. (pq)-(|39f), we find after some algebra 
GJ cosi ,„ „ 2 . 2 



c 2 p 3 n 



2 tan 2 i cos 2 <p ). (40) 



-[li-y^-— (1-2 tan 2 icos 2 <p )], (41) 



As expected, e is proportional to J and vanishes for a nonrotating source. By means of 
eq. (|34]) it finally follows that 

2tt GJ cosi 

n c^p^ n 

or 

T + -T_ = 47T^^(1 - 2 tan 2 * cos 2 <p )- (42) 

c 2 M 

Let us note that for i = 0, we recover the result of Section 3. On the other 
hand, e diverges for % = ir/2. It follows from general relativity that for a geodesic 
(spherical) polar orbit, the clock effect disappears since the angular momentum vector 
of the source in effect lies in the orbital plane. The orbital period is then simply given 
by the gravitoelectric Keplerian period. On the other hand, the period for azimuthal 
closure is given by 27r(2GJ/c 2 p 3 ) -1 , which is very long compared to the Keplerian period 
fl3| . This circumstance is reflected in our first-order perturbative result given by eq. 
(41): For i — > 7r/2 , T± — > oo. Therefore, in eq. (41), the inclination angle i must 
be sufficiently less than tt/2 such that the perturbative treatment in Section 4 remains 
valid. 

An important feature of eq. (41) is that when i ^ 0,T± depends upon <p , i-e. the 
clock effect depends in general on the position of the mass m along the orbit at t — 0. 
This dependence of the clock effect on where the mass m is along the orbit when the 
timing observations begin is illustrated in Fig. 2 and could be helpful in the detection 
of this effect. 

Up to now we have assumed that when unperturbed the two satellites orbit along 
opposite directions in the same plane with arbitrary inclination i. Let us now consider 
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the case of two masses m + and m_ having the same distance from the center but moving 
in different orbital planes, say < i + < tt/2 and tt/2 < i_ < tt, respectively. Following 
the same reasoning as before, their orbital periods, as viewed by a static, asymptotically 
inertial observer and after choosing suitably t so that ipo = tt/2, can be written as 

T _= T (o) +2 ^_ cosi (43) 

c z M 

r=T (o)_ 27r _L cos ^ (44) 

c z M 



From these expressions, we immediately find 

J 



T_ -T_ = 2tt-— (cosz+ + cosz_), (45) 



T_+T_ = 2T (0) + 2tt— - -(cosi + - cosi_). (46) 

c z M 

From eq. (^) it is seen that the sum of the orbital periods of the two point 
masses will also show a gravitomagnetic contribution provided the inclinations of the 
two satellites are different. Note that eq. (^) reduces correctly to eq. (|42| ) if the two 
orbital planes coincide, while in eq. (pEBj ) the gravitomagnetic contribution vanishes. 
Moreover, a very interesting feature arises for supplementary inclinations of the two 
satellites, i.e. i + + z_ = 180°; indeed, in this case eqs. (^3|)-(^B|) become 



T_-T_= 0, (47) 

T_+T_ = 2T (0) + 4tt-^-cosz + . (48) 

c 2 M 



One of the most striking implications of eq. (|4*8|) is that the LARES mission ['I ! 
originally proposed to detect the Lense-Thirring drag of the orbital plane, could also 
be used to detect the gravitomagnetic contribution to the sum of the orbital periods. 



However, due to the present uncertainty of the value of GM ffi , <5(GM ffi ) = 8x 10 cm 3 s 



11 OTTT.3 „-2 



15|, the error in the unperturbed Keplerian period of the LAGEOS satellite is larger 



than the effect to be measured 

5T {0) = 1.52 x 1(T 5 s, (49) 

4tt^— cosi + = 4.71 x 1(T 8 s. (50) 

c 2 M 

Therefore, an improvement of our knowledge of GM^ would be necessary before 
any observables involving the sum of the unperturbed orbital periods and of the 
gravitomagnetic corrections may become detectable. 

Finally, these results suggest that the gravitomagnetic clock effect may be enhanced 
by considering suitable constellations of satellites orbiting the Earth. 
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5. Conclusions 

Exploiting the formal analogy between the law of motion of a charged particle acted 
upon by an electromagnetic field and the weak-field and slow-motion approximation of 
general relativity ("gravitoelectromagnetism"), it has been possible to derive in a simple 
fashion the gravitomagnetic clock effect for a couple of point masses following spherical 
orbits in space. General expressions for arbitrary values of inclination angle to the 
equatorial plane are given for these spherical orbits. These results have been obtained 
by neglecting terms of 0(c~ A ); in this way, eq. fl42|) is equally valid for proper periods 
of comoving clocks. If the two satellites orbit in planes with different inclinations, the 
sum of their orbital periods also exhibits a gravitomagnetic part. 

It is worth-while to compare our main result eq. ([4^) with the corresponding general 
relativistic expression (see p. 143, eqs. (36)-(37)), 

T i — T_ = 47r - - A cos a, (51) 

+ C 2 M ' V I 

where 

A' = A - 3$ r = T - 2Tq 1 tan 2 a cos 2 r] - 3$ r . (52) 

Here, a and r/ must be replaced by i and <p , respectively, and r — > 1, $ — » 
provided that terms of 0(c~ A ) are neglected. For the proper periods of comoving clocks 
an equation similar to eq. ( [51]) holds, except that A must be replaced by A. Note that 
these results have been obtained in a perturbative way and hold only for a sufficiently 
different from it/2. Thus our approach gives the same result for T + — T_ as general 
relativity once terms of order c~ 4 are neglected. 

As it is for all general relativistic effects, the observation of the gravitomagnetic 
clock effect is a tremendously difficult undertaking as well. This becomes immediately 
clear by noting that for a near-Earth orbit (T ~ 10 4 s) a time variation of ~ 100 ns can 
be caused equally well by a radial or azimuthal deviation of ~ 0.1 mm from the ideal 
orbit and therefore all forces that may produce accelerations larger than ~ 10~ 12 m/s 2 
must be taken into account. The empirical verification of the clock effect essentially 
faces two problems: (a) to measure with the utmost precision the actual position of the 
satellites, (b) to model with extreme accuracy all perturbing forces which will influence 
the period of the satellites. For a single orbit, this goal is certainly unattainable, but it 
may become feasible after a sufficiently long time of observation due to the accumulative 
character of the clock effect. While the present satellite-to-satellite tracking techniques 
allow the determination of an orbit with an accuracy of ~ 1 cm so that a minimum 
of ~ 1000 revolutions will be needed for the clock effect to become detectable, the 
consideration of the perturbing forces at the required level is quite demanding. Non- 
gravitational perturbations may be overcome by means of modern drag-free technology; 
however, the correct determination of all gravitational effects is limited by the accuracy 
of the respective Earth gravity field models that are currently available. Preliminary 
results suggest that it is in particular the uncertainty in the even zonal harmonics of 
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the spherical expansion of the terrestrial gravitational field as well as the zonal 
tidal perturbations [|| that presently inhibits the successful realization of the clock 
experiment. It should be mentioned, however, that upcoming geodetic space missions 
(especially GRACE and GOCE) are expected to improve significantly the accuracy of 
the gravity field of the Earth and may then allow the observation of the gravitomagnetic 
clock effect within a few percent accuracy 
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Figure 1. Schematic diagram of the unperturbed circular trajectory for the general 
spherical orbit. In practice, the instantaneous orbital plane at the instant observations 
begin, i.e. t = 0, can be taken to be the plane of the unperturbed orbit. 
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Figure 2. Plot of AT = T + — T_, given in eq. (42), in seconds versus the orbital 
inclination i in degrees for a satellite in a spherical orbit around the Earth. For inclined 
orbits, the three graphs illustrate the dependence of the gravitomagnetic clock effect 
on cpo, which is the angular position of the satellite along its circular orbit at t = 
measured from the line of the ascending node (cf. Figure 1). 



